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Abstract. At the light of recent results in literature we review a con- 
jecture formulated in Math. Phys. Electron. J. 1 (1995), paper 5, 1-13, 
about the mechanism of breakdown of invariant sets in KAM problems and 
the identification of the dominant terms in the perturbative expansion of 
the conjugating function. We show that some arguments developed therein 
can be carried out further only in some particular directions, so limiting a 
possible future research program, and that the mechanism of break down of 
invariant tori has to be more complicated than as conjectured in the quoted 
paper. 



1. Introduction 



1.1. The state of the art of [GGM]. In [GGM] a conjecture about the mechanism of breakdown 

of KAM invariant tori is proposed. Roughly it is based on the following idea (we refer to [GGM] 
for a more detailed and technical exposition and also for the introduction of the notions used 
in the following analysis). 
If T = R/27rZ, let be the ^-dimensional torus. Gonsider a Hamiltonian system 

n=ujj^- A + ^A- J-^A + ef{a) , (1.1) 

where (a, ^1) € T x M are conjugate variables, J is the matrix of momenta of inertia, • denotes 
the inner product in R^, /(a) is a trigonometric polynomial in the angle variables and e is a 
parameter. 

For concreteness one can suppose that TL describes the Escande-Doveil pendulum, [ED] : £ = 2 
and /(a) = acosai + 6cos(ai — 02), with (a, h) e R^. 

The solutions of the equations of motion describing the invariant tori for the system (1.1) 
with Diophantine rotation vector a; can be parameterized as 

a = ±+h{±,e), A = Aa + H{ir,e), (1.2) 

where ij^ G T'^, Aq = J{uj_ — Wg) and /i, i?, for e small enough, are analytic functions, whose 
series expansions admit a graph representation in terms of trees; we refer to [GGM] for details 
and definitions. 

One can consider also trees in which no resonances (see [GGM], §4) are allowed to appear 
but the perturbative parameter e is replaced by = e(l — o-g)"^, with a suitable (matrix) 
form factor taking into account the resummation of all resonances, [GM], and denote by 
h*,H_* the functions so obtained; the series expansion in 77 will be called resummed series. 

If p is the radius of convergence of the series defining in (1.2) and £c S is the 

(positive) critical value at which the tori break down, one has £c > P'-, in general the analyticity 
domain in e of the scries defining h, H_ is not a circle and it can happen that Ec > p- 

One can imagine that the following scenario arises: the singular behaviour of h{tp; e),H_{tl!; e) 
as functions of s is the same of that of h*{ip;r]e),H_*{tp;r]e) and while, for e — > e~, is still 
finite, r] = rj^ goes out the convergence domain (in t/) of the series for h*,H_*. If moreover the 
analyticity domain in 77 of /i* , H_* turns out to be a circle this would mean that 77 is a more 
natural parameter for the invariant tori. 

If this really happens then the behaviour of the series near the critical value is determined 
by the trees without resonances, for which the perturbative series would still be meaningful for 
£ near Sc- 

Under such an assumption a universal behaviour of the series near the critical value is pro- 
posed in [GGM]. Predictions of the value of the critical exponentc 5 are made, conjecturing that 
only some simple classes of trees are relevant: the linear trees presenting the largest possible 
number of small divisors (see [GGM], §5). 

Considering (for concreteness purposes) a rotation vector w = (r, 1), where r is the golden 
mean (sec §3.1 below), then one was led to expect that, by denoting by the denominators of 
the convergents defined by the continuous fraction expansion for r, a suitable function Z{Aq^), 
defined in equation (5.2) of [GGM] and recalled in §3.1 below, satisfied the asymptotics 

C = C{r„f), (1.3) 
for some positive constant C{r], f) depending on r]f. 

1.2. About the standard and semistandard maps. If true, the above mechanism should work 
also for the standard map, which is the dynamical system generated by the iteration of the 
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area-preserving map of the cyHnder to itself 

x' = X + y + esinx , (14) 
y' = y + esinx , 

where (x, y) G T X M. At least it is generally assumed that arca-prcscrving maps and Hamilto- 
nian flows share the same critical behaviour: in particular they are expected to have the same 
critical exponent 6; see also [GGM], [CGJ] and comments therein about [Ml], [M2]. 

The role of the functions h, H_ is now played by two scalar functions u, v (see, for instance, 
[BGl]) such that 

x = a + u{a,s), y = 2Truj + v{a, e) , (1-5) 

with u) a fixed rotation number and u, v analytic in their arguments ^. In terms of a G T, the 
dynamics is a trivial rotation: a — > a' = a + 2Trw. 

Analogously to the Hamiltonian case (1.1), one can define two functions u*,v*, obtained by 
considering only the values of the trees contributing to u, v without resonances (and replacing 
£ with 7], where r] = % = s{l — (Te)~^ takes into account the resummation of resonances; note 
that now is a scalar). 

For the standard map, let us denote by p{uj) and ec(t^), respectively, the radius of convergence 
(in e) and the critical value for fixed rotation number lu. 

Consider also the semistandard map, introduced by Chirikov, [C], 

( x' = X + y + e' expix , . . 

\y' = y + exp ix , s' = e/2 , \ ■ ) 

and call uo,vo the corresponding conjugating functions; let us denote by po(w) the radius of 
convergence (in e') for the semistandard map^. 

1.3. The implications of [Dl]. Given a rotation number co G (0, 1), let B{u)) be the function 

Bi.)=j:'^, (1.7) 

n=0 ^" 

where {qn} are the convergents defined by the continued fraction expansion for the rotation 
number u = [ai, a2, a^, ■ ■ ■], so that (j„ = a„g„_i + (j„-2, for n > I, and g_i = 0, Qq = 1. The 
function (1.7) is related to the Bryuno function introduced by Yoccoz, [Y]. 
Consider the series defining u* = u*{a, rj). Of course one can write 

oo / k \ oo 

u*{a,v) = Y^n" E = Er?VW(a) . (1.8) 

fe=l \v=-k ) k=l 

The radius of convergence R{co) of a series like (1.8) can be expressed as 



ae[0,27r] fe— >oo 



1/fe 



R ^(w) = sup limsup M*^''' (a) = limsup max u*^''^ , (1.9) 



fe^oo kl<fe 



i/fe 



as the equivalence between the two definitions has been shown in [Dl]. 

On the other hand the coefficient in (1.8) with i/ = k is the same for both the functions u 
and u* (and one has u^,^^ = u*^^'' = Uq'^). Moreover it is easy to see that one has 



,(fe) 



> B^e'^^^'^^'' (1.10) 



^ The functions u,v axe trivially related: v{a,e)=u{a,e)—u{a—2wu!,£). 
2 

We could define also the functions u^.v^ and the critical value eoe('^) for the semistandard map, but one has 
uo=Ug,vo=Vg, as the semistandard map has no resonances, and eoc(.'^)=Po{'^)- to deduce the latter property 
simply note that the semistandard map is invariant under rotation of e' in the complex plane. 
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i; again see [Dl]. 



*(fc) 



i/fe 




(1.11) 



so that R{ij) < B^^e 2-8(1^)^ j g ^j^g radius of convergence in r] of the function u* at best 
should be of order of po('^) ^• 

1.4. The implications of [D2]. One has po{io) = Ci(ti;)e~^^^'^) with Ci(a;) satisfying the bound 
Cj~^ < Ci{lo) < Ci, uniformly in cu, for a suitable constant Ci, [Dl]. Also for the standard 
map one has the same dependence of the radius of convergence p{u)) on the function B{u), with 
a possibly different function Ci{co), always admitting a bound from below and from above: the 
bound from above is proven in [Dl], by using the argument recalled in §1.3, and the one from 
below in [BG2]. 

In [D2] it is proven that for rotation numbers ^ w = 7„ = [n, n, n, . . .], with n large enough, 
one has C2/n > £c(7n) > /n, for some constant C2; as B{u)) > logn, for u) = 7„, then 
^dln) > C2^n~^ and p{^n) < Cin~^, (see [D2], end of section 5). Therefore one can choose 
n so large that p{jn)/£cijn) is smaller than any prefixed quantity. 

Note that the difference between the radius of convergence and the critical value becomes 
relevant only for some rotation numbers (like the noble numbers 7„, with n large enough). 
For instance for 71, according to numerical simulations, the analyticity domain appears (very) 
slightly stretched along the imaginary axis so that not only the two quantities are comparable, 
but even £c(7i) = P(7i), [FL]. 



2.1. About the (1-3). If the dominant contributions to were given by summing only the 
values of trees defining the functions Z(Ag^), as conjectured in [GGM], then the standard map 
and the semistandard map should have the same critical behaviour. In fact in terms of trees 
the semistandard map admits the same graph representation of the standard map, with the 
only (remarkable) difference that the mode labels have all the same signs; then the paths with 
the maximal number of small divisors in Z{Kq^ ) are the same for both the standard map and 
the semistandard map. 

This means that, accepting the above conjecture, one ought to expect that the radius of 
convergence (in e) of the functions uo,vo for the semistandard map and the radius of conver- 
gence (in T}) of the functions u* , v* for the standard map should be equal to each other and, in 
particular, of size of Pq{u}). 

Furthermore the results listed in §1.3 imply that the functions u*,v* can not have a radius 
of convergence larger than that of the functions uo,vo, also considering all trees contributing 
to them and not only the ones defining the function Z(Ag^). 

In other words the semistandard map should capture all the critical behaviour of the standard 
map: if so there would be contradiction with the results existing in literature. As a matter of 
fact we shall see in §3 that the behaviour (1.3) conjectured for Z(Ag^) does not hold. 

2.2. About the behaviour of near Sc- The idea that the critical behaviour can be studied 
through the series obtained be neglecting the resonances (that is by resumming them and 
defining a new parameter 77 = rje, as briefly recalled in §1.1 and in §1.2) is presented in [GGM]. 
One could also interpret the numerical results of [CGJ], [CGJK] as supporting such an idea: 
the fact that the KAM iteration represents the good transformation to look at also far from 

o 

In the same way one finds that the radius of convergence in e of the function u at best should be of order of 
^ With the notations in §1.1, one has r=7i. 



2. Discussion 
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the KAM analyticity domain could suggest that a perturbative approach to the study of the 
breakdown of KAM invariant tori in some sense should be possible: this was the idea underlying 
the analysis performed in [GGM] . Anyway much work has still to be done in this direction; see 
also §4 below. 

From the results listed in §1 the following picture emerges. For e small enough < R,\e\ for 
some constant R; see the theorem in [GGM], §4. When e grows along the real axis toward the 
critical value, the possibility to have still finite for e = Sc is consistent with a perturbative 

approach if one of the following two cases arise: 

(1) if cTg is finite and smooth in e for e = £c, then it can happen that jryd = |£c(l — '^So)~^\ 
becomes equal to the radius of convergence R{u}) of the functions u*,v*; 

(2) if CTe is finite and singular in e for £ = Sc, while rjc = £c(l — Cec)"^ is smaller than the 
radius of convergence R{u!) of the functions u*,v*, then the singularity of the tori shows up 
through the singular dependence of de on s at the critical value. 

In both cases the perturbative series for u*,v* can be used also near the critical value, because 
for any |£| < Sc one has that rj is smaller than the radius of convergence of the series. 
In the case of a; = 7„, in principle two possibilities can be envisaged for the behaviour of ae 

for £ near Sc- 

(i) cither \(t^\ < 1, 

(ii) or else negative and \ae\ > 1 ^. 

Consider lu = jn, with n large enough, and assume the results in [D2]. Then, if the analyticity 
domain in r] was a circle with radius R{uj) and, for s s~{oj), \ae\ < 1 and rj -R(w), one 
should have |?7| = |£(1 — <Te)~^\ — > R{u)) « po{u)), hence |1 — ae\ — > |1 — a^J ~ £c(w)/po('*^) > 
Con, with Cq^ = C1C2, which is incompatible with |cr£| < 1. So only the possibility (ii) above 
could be consistent both with [D2] and with the scenario proposed in [GGM], for w = 7„. 

More generally, for u = 7„, with n large enough, if one had r] ^ rjc < R{io) for e e7(w), 
then \ae\ should become at least of order e^^'^^ for s s~{uj), i.e. \ae\ \aej > C3{u!)e^^'^\ 
for some function Cs^uj) bounded uniformily in lo: in fact in this way one would have [r^d = 
|£c(l — Ue^)"^! ~ C3£c(w)e~^('^) C4po(a;) = -R(w), for some constants CsjC^, in case (1) and 
1 77c I < C4Po{lo) = R{u)) in case (2). 



3.1. The definition of Z{Aq^). If r = (\/5 — l)/2 is the golden mean, let us call {pn/Qn} the 
convergents of the continuous fraction expansion for r, where {p„} is the Fibonacci sequence 
defined by Pn+i = Pn + Pn-i with p-i = 1 and po = 0, so that pn = Qn+i with g_i = and 
qo = 1- 

The sequence of numbers Z„ = Z{Aq^) is defined in [GGM], equation (5.2), as sum of the 
values of a suitable class of trees which can be described by the family Ag^ of self-avoiding 
walks on Z starting at (0,0), ending at {qn,—Pn) and contained in the strip < x < qn, 
except for the left extreme points. 

Then the numbers Z„ can be approximately defined by the recursive relation 



where, for consistency, we fix Z_i = and Zq = r ^, [GGM]. Set also £„ = qnV — p„; then 
£-1 = -1 and £„ = g„r - g„_i. 

^ Note that a., positive and |<Te|>l is not possible as the analyticity of ae in e for small e would imply the 
existence of a value e such that |1— <75|=0. 



3. Analytic results 




(3.1) 
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Instead of studying (3.1) define A„ = log Zn, so that 

A„+i = A„ + Xn-l + 2 (log£„_i - log£„+i) , (3.2) 

with A_i = 2 log r and Ao = — 21ogr. 

3.2. Against the (1.3). Ah the following identities can be easily proven by induction using the 
fact that r is the positive solution of + r — 1 = 0. 
First note that e„ = (-l)"r"+^ It follows that 

A„+i = A„ + A„_i-41ogr . (3.3) 

which satisfies A„ = Atn + Snj where 



,n — 1 

So = s_i = , 



s„ = -41ogrX;i=o ' forn>l, ^3^^ 



and fjtn is defined by = /x„ + Hn-i, with = 21ogr and /xq = — 21ogr. It is now 

immediate that Hn = — (21ogr) qn-2, for n > 1, so that we have 



n-l 



Xn = -2 logr Un-2 + 2 ^ (?, , (3.5) 



i=0 



for n> 1. 
Let now use that 

^-(n+l) _ ^_]^^n+lj.n+l 



Qn = • (3.6) 

r + r ^ 



Inserting this expression in (3.4) we get, for n > 1, 



A„ = ii^^ (r-("-i)+2yr-'-(-l)"-V"-i-2y(-l)vl = 

^ + ^-' V h h J (3.7) 

= 21ogr-i (r-("+2) _ 2 + (-l)"+2r"+2) , 
which can be written as 

Xn = 21ogr-i [(r + r-^)qn+i - 2 + 2(-l)"+V+2] , (3.8) 

which means 

Zn = r\-2ir+r-')<ln+ir-^i-ir+'r"+- . (3.9) 

This implies, using that 



j.n+1 ^ 1 -1 



l + (-l)"r2("+i)l , (3.10) 



the following asympotics: 

\Z{A,J\ « KC-^^ (1 + i-irdq-l,) , (3.11) 

with K = r~^, C = r~^^^'^^ and d = 4(r + r^^)^^ logr~^. This is clearly in contradiction 
with (1.3). 



4. Conclusions 

4.1. Dominant contributions. Our attention to [GGM] has been called back by the recent pa- 
pers [CGJ], [CGJK], where the breakdown of KAM invariant tori (for two-dimensional Hamil- 
tonian systems) is numerically studied through a renormalization group scheme, and by the 
results in [D2]. 
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Even if it is not true that the terms defining Z{Aq^) are the only relevant ones (as the results 
in §3 show), one can argue that the real dominant contributions are given by the trees having 
the mode labels accumulating near the resonant line {i.e. such that the small divisors are really 
as small as possible), but not necessarily belonging to the class described by A^^. 

In fact let us compare the tree values for the Escandc-Dovcil pendulum with the ones for the 
standard map and accept that they have the same singular behaviour at the critical value (as 
it is generally believed, [ED]): the small divisors are respectively (iw • v)'^, v = (z^i, ^2) € Z^, 
and 2[cos(27ra;i/) — 1], v G Z. So one can note that they have the same smallness problem, 
with the only difference that for the Escande-Doveil pendulum one can have also small divisors 
which are not small at all (when u is nearly parallel to of), while for the standard map one 
can approximate the quantity 2[cos(27ra;z^) — 1] with something of the form {i{wi'i + 1^2))^, 
with i^i = ly and 1^2 such that {uji'i + 1^2! < 1: in other words one can expect that the tree 
expansion for the standard map is very similar to that of the Escande-Doveil pendulum, but 
it does contain only the trees with the momenta directed along the resonant line. Of course 
this imply neither that in the case of linear trees the most dominant contributions arc the ones 
with the mode labels having all the same signs nor that only the linear trees are relevant. As 
a matter of fact the analysis performed in §3 shows that making a so restrictive assumption 
leads to wrong results. 

4.2. Resummed series: smooth form factors. In conclusion an overall behaviour like (1.3) 
could be still possible in principle, even if a larger class of trees ought to be taken into account. 

Let us consider oj = 7„, with n large enough and use the results in §2.2. If for e ^ Sc one had 
|ae I < 1, then we have seen that the analyticity domain in t] oiu*, v* can not be a circle: rather, 
defining r]c{u>) as the equivalent of Sc for the functions u*,v*, the ratio r]c{i^) / R{u)) would be 
of the same order of the ratio ec{<-^)/ pii-^), and no simplification could arise in considering the 
resummed series u*,v* instead of the original ones u,v. In other words rj would not be a 
"natural" parameter. 

Of course in the case of the standard map, deep cancellation mechanisms should intervene 
in order to enlarge the analyticity domain along the real axis of the series for u* , v* from 
i?(w) « po{u)) to a quantity rjdco) of the same order of ec(i^)- Analogous considerations could 
be made for the functions ht ,H* in the case of Hamiltoiiian flows. 

On the contrary if, for e Ec, one had case (1), possibilty (ii) - sec §2.2 so that 
C2(u;)e^("\ the convergence domain in rj for u* , v* could still be a circle. Of course, if this is 
the case, the study of the factor form could turn out to be a very difficult task: in fact it 
could become a nonperturbative problem, as the full dependence on e would be required for 

£ Sc. 

A perturbative analysis could still be possible for some rotation number, for instance for 
w = 7i, when the radius of convergence and the critical value are expected to be equal (see 
comments at the end of §1.4); anyway we have no evidence of this, so we leave it as an open 
problem. 

4.3. Resummed series: singular form factors. Suppose now that at the critical value is finite 
and singular. One can consider the resummed series for the form factor CTe, by expressing it as 
series in 77, i.e. = F{rj) = F(£(l— Ue)"-^), where the function F admits a graph representation 
in terms of trees without resonances and with e replaced with 77^. If for e = Sc one has that rj is 
inside (enough) the domain in which the perturbative series for the function F converges, then 
one can try to truncate the series to the first orders, so obtaining an (approximate) implicit 
equation for a = a^: the solution should be singular at the value e = Ec- So the possibility of 
a perturbative study of the breakdown of KAM invariant curves has not to be excluded ^. 

Further investigation (also from a numerical point of view) in this direction would be highly 

^ Note that for w=7„, with n large enough, one can still define the function F{r]): the problem is that for e 
near to the critical value no truncation of the series would be meaningful. 



6 



profitable. Also a comparison between the values of the involved quantities for the standard 
map and the semistandard map would be very enlightening. 

Acknowledgments. We want to thank the ESI in Vienna, where part of this work has 
been done, for hospitality. Wo also thank V. Mastropietro for interesting discussions and in 
particular G. Gallavotti for continuous enlightening suggestions and critical comments. 

Appendix Al. Numerical analysis 

Al.l. Motivations. Although the results of §3 are conclusives we report here on some numerical 
simulations that suggested us those results. We tried to fit A„ with a slightly more general 
relation than (1.3), i.e. A„ = + 5 logg„ + cg„. The problem of this fit is that g„ ~ r~" ,so that 
we have to compute all the involved quantities with a very high accuracy: if not the constant 
k and the "linear" term S will be completely lost. We think that the way we did can be of 
some interest for the reader. 

We used the standard Unix command be that is able to execute computation written in a 
C-like programming language that operates in fixed point notation with number of arbitrary 
dimension and an arbitrary but prefixed precision. We compute Cj, fcj, 6i has the best square fit 
of (3.1) using the value of Xm with m= (i — 1)100 to m = ilOO. The only problem is to choose 
the precision at which the computation are done in such a way to have a desired precision in 
the value of Ci,ki,Si. This is what we will discuss in the next subsection. 

A1.2. About the preeision. Let p be our chosen precision, i.e. we make all operation with 
p significant digits after the point. This mean that at every elementary operation creates an 
error 0(10"^). 

If we call 9i the error due to this round off at step i it's easy to see that the accumulated 
error on A„ is X]r=o (^iQn-i-i — O(10^Pr"). To compute the best quadratic fit we have to 
solve the equation Ax = v, where A is the symmetric matrix A = X^"^o *5n ^ Qn, with 
Qn = (g„,log((7„),l), V = YTi=aQn>^n and x = {c,S,k). This mean that Ai^i « r"^", while 
Ai^2, Ai^3 « r~" and all the other entries Aij of A, with i > j, are of order 1. 

Observing that {A-'^)ij = A^'-'/det A and that det A w r'^" ± IQ-^r-^", we get 

' (A-i)i,i « r^" ± IQ-P , 

< (A-i)i,3,(A-i)i,2~r"±10-fr-", (4.1) 
^ (A-i)2,2 , (A-i)2,3 , (A-i)3,3 ~ 1 ± lO-fr-2" . 
The above estimates with the fact that u w (r-^^ilQ-Pr-^", r-^ilO-Pr"", r-^ilO-Pr"") 
implies an error on x of order 10~^r~^"'. 

A1.3. Num.erical results. The numerical result that we get are summarized in the next table, 
which, in our opinion, requires no comments. According to our error analysis, we fixed the 
precision to p = 2000 log^g + 20 in such a way that the reported numbers should be reliable 
within 20 digits. The agreement with the analytical evaluation reported in §2 is perfect. 



n 


c 


5 


k 


100 


3.482081477708027334 


-.000312407255699678 


-1.915302823598781109 


200 


3.482081477708027334 


-.000000000000000000 


-1.924847300238413790 


300 


3.482081477708027334 


-.000000000000000000 


-1.924847300238413790 


400 


3.482081477708027334 


-.000000000000000000 


-1.924847300238413790 


500 


3.482081477708027334 


-.000000000000000000 


-1.924847300238413790 


600 


3.482081477708027334 


-.000000000000000000 


-1.924847300238413790 


700 


3.482081477708027334 


-.000000000000000000 


-1.924847300238413790 


800 


3.482081477708027334 


.000000000000000000 


-1.924847300238413790 


900 


3.482081477708027334 


,000000000000000000 


-1.924847300238413790 


1000 


3.482081477708027334 


.oooooooooooooooooo 


-1.924S4730U238413790 
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Table 1. Numerical results. The small deviations for i=l is due to the fact that 
our fit for A„ contains g„ instead of qn+i- 
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